Thermoelectric Behaviour Near Magnetic Quantum Critical Point 
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We use the coupled 2d-spin-3d-fermion model proposed by Rosch et. al. (Phys. Rev. Lett. 
79, 159 (1997)) to study the thermoelectric behaviour of a heavy fermion compound when it is 
close to an antiferromagnetic quantum critical point. When the low energy spin fluctuations are 
quasi two dimensional, as has been observed in YbRh2Si2 and CeCue-a: Au^;, with a typical 2d 
ordering wavevector and 3d Fermi surface, the "hot" regions on the Fermi surface have a finite 
area. Due to enhanced scattering with the nearly critical spin fluctuations, the electrons in the hot 
region are strongly renormalized. We argue that there is an intermediate energy scale where the 
qualitative aspects of the renormalized hot electrons are captured by a weak-coupling perturbative 
calculation. Our examination of the electron self energy shows that the entropy carried by the hot 
electrons is larger than usual. This accounts for the anomalous logarithmic temperature dependence 
of speciflc heat observed in these materials. We show that the same mechanism produces logarithmic 
temperature dependence in thermopower. This has been observed in CeCue-aiAui,. We expect to 
see the same behaviour from future experiments on YbRh2Si2. 



I. INTRODUCTION 

Understanding the behaviour of a system close to anti- 
ferromagnetic quantum critical point (QCP) is currently 
an area of active research. The problem is interesting 
both in the context of high temperature superconductors 
as well as heavy-fermion materials, especially to under- 
stand metallic phases that show non-Fermi liquid (NFL) 
properties. In recent times several materials have been 
discovered where it has been possible to demonstrate the 
existence of magnetic QCP.^"'^ This has made the prob- 
lem an exciting ground where theoretical understanding 
of electrons with strong correlations can be verified ex- 
perimentally. One central issue in this problem is an ap- 
propriate theoretical treatment of electrons interacting 
with spin fluctuations close to the QCP where magnetic 
correlation length diverges. A second central issue, is 
whether the spin fermion model describes the relevant 
degrees of freedom, or whether a more basic model, al- 
lowing for the disintegration of the binding of local mo- 
ments to the quasiparticles, is necessary for describing 
this transition. ^'^ 

In this paper we will discuss two experimentally 
well-studied heavy fermion materials, CeCug-xAua:^ 
and YbRh2Si2,^ that exhibit antiferromagnetic QCP. 
In doped CeCug, replacing Cu with larger Au atoms, 
favours the formation of long range magnetic order. ^ Be- 
yond a critical doping Xc = 0.1, the ground state of the 
system is antiferromagnetic with finite Neel temperature 
(Tjv).^ At the critical doping Tjv is zero and the system 
has a QCP. On the other hand YbRh2Si2 is undoped 
and atomically well-ordered.^ It is a much cleaner ma- 
terial than CeCvLe-xAux, with residual resistivity (po) 
smaller by a factor of about 10. At ambient pressure it 
develops long range magnetic order at a very low tem- 
perature of Tpf ~ 65 mK.^ The ordering temperature can 
be suppressed to practically zero (less than 20 mK) by 
applying a magnetic field of only 45 mT.^ Both these 



materials show pronounced deviations from Fermi liquid 
(FL) behaviour, which is beheved to be due to closeness 
to the QCP. For instance, the dependence of electrical re- 
sistivity Ap = p — pq to temperature T is Ap oc T, while 
that of specific heat C is C/T oc -InT.^'^ This is in 
contrast with FL behaviour which predicts Ap oc and 
C/T = constant. The low temperature NFL behaviour 
is observed over a decade of temperature, up to about 
1 K for CeCug-a^Aua;,^'^ and up to as high as 10 K for 
YbRh2Si2.^ The source of the interesting physics in these 
materials is the localized 4/ electrons^ of Ce'^+ (in 4/^ 
electronic configuration) and Yb'^^ (in the configuration 
4/^^), and their interaction with the relatively delocal- 
ized s, p and d orbital electrons that form a conduction 
band with a well defined Fermi surface at low tempera- 
ture. The conduction electrons and the localized 4/ elec- 
trons carrying magnetic moment are coupled by exchange 
interaction (J) . Below a certain critical value of exchange 
interaction (Jc), the local moments interact with each 
other, mediated by conduction electrons, and at suffi- 
ciently low temperature form long range antiferromag- 
netic order. On the other hand, if the exchange coupling 
is strong (J > Jc), the local moments are quenched below 
a certain temperature (lattice Kondo temperature) . The 
quenched moments hybridize with the conduction elec- 
trons and they participate in the formation of the Fermi 
sea. The ground state of such a system is non-magnetic. 
The exchange coupling is usually tuned experimentally 
by either doping the material or by applying external 
pressure or external magnetic field. 

For CeCug-xAux there are two different views'* regard- 
ing the nature of the system in the non-ordered phase and 
the corresponding mechanism by which the critical insta- 
bility occurs. In the first picture, the lattice Kondo tem- 
perature (T^) becomes zero exactly at the critical point 
(J = Jc). The local moments of the 4/ electrons survive 
at all finite temperature close to the critical point. At the 
transition point they are critically quenched. The local 
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moments produce the critical magnetic fluctuations that 
destabilize the Fermi sea. It has been argued, in favour 
of this mechanism, that the data on magnetic suscep- 
tibility shows non-trivial scaling with temperature/ At 
the critical point the susceptibility has the scaling form 
^ _ T~"f{uj/T) with an anomalous exponent a ~ 0.75, 
which is different from conventional insulating magnets 
which have a= 1. The alternative picture suggests that 
is finite at the critical point. Well below this tem- 
perature, and close to the critical point, the local mo- 
ments are quenched by Kondo mechanism. The 4/ elec- 
trons become part of the Fermi sea. Then, the phase 
transition occurs by the usual spin-density wave insta- 
bility of the Fermi surface. In this picture the local mo- 
ments do not play any role in the phase transition. This 
theoretical viewpoint, proposed by Rosch and collabora- 
tors.^ is motivated by inelastic neutron scattering data 
on CeCus.gAuo.i which show that the nearly critical spin 
fluctuations are two dimensional.^ But the origin of the 
quasi 2d behaviour of spin fluctuations is not well un- 
derstood. However, the same feature is probably also 
present in YbRh2Si2, where the structure of the lattice 
provides a more natural explanation for the spin fluc- 
tuations to be 2d.^ Besides the nature of the magnetic 
correlations, there are different opinions regarding the 
dynamics of the spin fluctuations. It has been argued^ 
that if the ordering wave- vector spans different points of 
the Fermi surface, then the dynamics of the spin fluctua- 
tions is overdamped, with dynamic exponent z = 2. This 
model of spin fluctuations with d = 2 and z = 2, coupled 
with three-dimensional electrons, was used to explain the 
linearity of resistivity with temperature.® Following the 
method of Hertz, in which the system is described 
entirely in terms of the spin fluctuations, after a for- 
mal Hubbard-Stratonovich transformation to integrate 
the fermion modes, it also explains the logarithmic tem- 
perature dependence of speciflc heat.®'^^ In an alternative 
description,^^ in terms of low energy electrons interact- 
ing with spin fluctuations, it has been suggested recently 
that both the frequency and momentum dependence of 
the spin fluctuation propagator undergo singular correc- 



tions such that the propagator acquires an anomalous 
dimension rj ^ 1/4.^^ Thus, after nearly a decade, there 
is still no clear understanding regarding the appropriate 
model that describes the quantum phase transition. 

In this paper we will study the thermoelectric be- 
haviour of a system in the paramagnetic phase and close 
to antiferromagnetic QCP. For CeCus.gAuo.i it is known 
that thermopowcr (St) has a dependence similar to spe- 
cific heat over the same range of temperature, ^^'^'"^ i.e, 
St/T (X —\nT. We will show that scattering with nearly 
critical spin fluctuations give rise to temperature depen- 
dent quasiparticle mass (m*) over much of the Fermi sur- 
face. The signature of this can be seen in static response 
(specific heat) and in transport (thermopower). Finally 
we will argue that the same mechanism should be rel- 
evant for YbRh2Si2, and so we expect to see the same 
behaviour for thermopower from future experiments. 



II. MODEL 

Our model is motivated by the second picture as de- 
scribed above. It assumes that defines a high energy 
parameter. For T ~ the local nature of the spins of 
the 4/ electrons is important as they participate in some 
lattice Kondo phenomenon. For T < T^, the 4/ elec- 
trons become part of the hybridized conduction band. 
In this regime the nearly critical spin fluctuations of the 
conduction electrons is important. It is an intermedi- 
ate temperature range where the system is described by 
low energy conduction electrons interacting with quasi 2d 
spin fluctuations. Within the spin- fermion description, 
at sufficiently low temperature, the 3d nature of the spin 
fiuctuations is retrieved and the model used here ceases 
to be valid. In this regime, the model predicts, in pure 
systems, a crossover to an electronic Fermi liquid with a 
finite mass. However the physics governing this dimen- 
sional crossover, has not been investigated. 

The model is described by the Hamiltonian 



n = X! ^fc'^U'^k.a + Y J2 Ck+q,aCk,/3a„,^ • S_q X [X ^(q)Sq-S_q-Fnq-n_q] 

k,CT k,q,c«,/3 q 

+T S (Sk,-SkJ(Sk3-SkJ^(ki + ... + k4). (1) 

ki,...,k4 
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Here cj^ ^ is the electron creation operator, Sq is the op- 
erator for the spin fiuctuations, Ilq = dtSq is the conju- 
gate momentum field for the spin fluctuations, and x(l) 
is the static magnetic susceptibility, go is the bare cou- 
pling between the electrons and the spin fiuctuations, and 
uo is the interaction energy of the spin fluctuations. The 



collective spin fluctuations are formally obtained by in- 
tegrating out high energy electrons in the band up to a 
certain cutoff. Thus the typical energies of the spin fluc- 
tuations LOs ~ W, the bandwidth of the conduction elec- 
trons. The system is close to an antiferromagnetic insta- 
bility with ordering wave- vector Q. We will assume that 
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the dynamics of the spin fluctuations is purely damped 
with dynamic exponent z = 2. The spectrum of the 2d 
spin fluctuations will be described by^°'^^ 



X 



"\q,u;) = 8 + Ws(q- Q) 



(2) 



Here b is the mass of the spin fluctuations and mea- 
sures the deviation from the QCP, the parallel directions 
are those along the planes of magnetic correlation, and 
7 ~ (ffo/cF)^ is an estimate of the damping from the 
polarization bubble. In the spin fluctuation part of the 
Hamiltonian, the interaction term Uq is marginal, since 
the scaling dimension is zero.^*^'^^ The main contribution 
of this term is to renormalize the mass of the spin fluctu- 
ations (5) and make it temperature dependent. Within a 
Gaussian approximation, b is linearly dependent on tem- 
perature, up to logarithmic corrections. ^'^^ We will ig- 
nore other effects of the term in our discussion, and 
will consider only the quadratic term with temperature 
dependent mass of the spin ffuctuations. To simplify the 
calculation we will assume a spherical Fermi surface for 
the non- interacting electrons, with the ordering wave- 
vector Q = (a, 0, 2fci? cos^^o)■ Here ^ (i.e. not 
Ikp ordering), and 7^ 7r/2 (i.e. not ferromagnetic or- 
dering). We have chosen x as the direction along which 
the spin fluctuations are uncorrelated, and a, the order- 
ing in the x-direction, varies from one plane of magnetic 
correlation to another. Since the spectrum of spin fluc- 
tuations is 2d, those carrying momentum of the form 
Q -h ax, where a is arbitrary, are all nearly critical. Due 
to constraints from energy-momentum conservation, only 
those points on the Fermi surface that arc connected by 
the nearly critical spin fluctuations are particularly sen- 
sitive to the QCP, since electrons at these points undergo 
singular scattering with the spin fluctuations. These are 
the so-called "hot spots". It is important to note that 
since the spin fluctuations are 2d, there will be a flnite 
area of the Fermi surface that is hot. Though it is worth- 
while to estimate the fraction of the Fermi surface that 
is hot, theoretically it is a daunting task. In our cal- 
culation we will assume that most of the Fermi surface 
is hot. In effect, we are assuming that contribution to 
static response and also to transport is mostly from the 
hot regions. It was pointed out by Hlubina and Rice^^ 
that in transport the hot carriers are less effective than 
the cold ones. This is because the quasiparticle lifetime 
of the hot carriers is less than that of the cold carriers, 
since the former suffer enhanced scattering with the spin 
fluctuations. As we will show below, the lifetime of the 
hot electrons cx 1/T, while the cold electrons have 
Fermi liquid characteristics with Tc l/T"^ . If x is the 
fraction of the Fermi surface (FS) that is hot, then we 
can make an estimate of conductivity cr. 

The first term, which is the contribution from the hot 



region, will dominate to give Ap oc T only if x > 
1/(1 -\-T /ep). This gives a rough estimate of the fraction 
necessary for the hot carriers to dominate. In the case of 
CeCue-xAua;, which is a dirtier material, the above es- 
timation is more involved. It was recently shown^'' that 
the effect of disorder is to favour isotropic scattering and 
thereby reduce the effectiveness of Hlubina-Rice mecha- 
nism. Thus, one should expect a smaller fraction, than 
estimated above, enough to make the contribution of the 
hot carriers significant for CeCug-xAuj,. 



III. ELECTRON SELF ENERGY 

To calculate the effect of the low energy spin fluctua- 
tions on the hot electrons, we will examine the electron 
self-energy. The lowest order term in perturbation gives. 



(p'-) = -|E 



2-Ki 



'-. x(k,f^)G(p + k,a; + 0), 



(3) 

where G (p, to) is the free electron propagator given by. 



G(p,a;) 



+ 



l-n„ 



CO — €p — IT] bJ — Cp + ir} 



Here Hp is the electron occupation of the momentum state 
p at T = 0. As expected, the above expression has dif- 
ferent behaviours in the hot and cold regions. But within 
each region the self-energy is practically momentum inde- 
pendent. The imaginary part of the self-energy gives the 
quasiparticle lifetime as determined by scattering with 
the spin fluctuations. For w > we have, 



ImE (p, (jj) 



V 



0<efc<a; U-|-Ws(k-Q)|| j -|- (w - Ck+p) 



If p is a point in the hot region, then it is connected to 
another hot spot by a wave- vector of the form k = Q-|-ax. 
We linearize the spectrum about this second hot point, 
and perform the integral in terms of local coordinates 
around it. In the hot region we get. 



ImS (p, w) oc 



( 9l 



yep^s) max[(5, ti)] 



(4) 



For u) > 5 the lifetime of the hot electrons is much smaller 
than that given by Fermi liquid behaviour (ImS(a;) oc 
uj"^). As we have mentioned above, this is due to more 
effective scattering with the spin fluctuations in this re- 
gion. For the cold electrons the behaviour is Fermi liquid 
Hke. 

Next, we will examine the real part of the self-energy. 
The dependence of ReS on frequency is more important 
than the dependence on momentum. We get. 
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lim ^ReS (p,w) = 



k+p 



(7^ 



In 



7k 



ek+p 



+ 



7r(2nk+p - l)7kek+p 



(.2 



■"k+p 



(5) 



Here 7k = + t<^'.s(k — Q)|. If p is a point within the hot 
region, each of the three terms in the above expression 
is logarithmic. As before, after linearizing the spectrum 
near the second hot spot, we get, 

- lim |-ReS (p, u;) oc ( In (^^) . (6) 



Due to scattering, the non-interacting electron mass m 
is renormalized to the quasiparticle mass m* = m/Z (in 
the absence of any momentum dependence of the electron 
self-energy), where 

d 

Z-^ = 1- lim — ReS (p,w) 

defines the quasiparticle residue. Since 5, which measures 
the deviation from the critical point, can be written as 
5 = T{p — pc) + T, the quasiparticle mass becomes tem- 
perature dependent. Here p is an experimental parame- 
ter that can be tuned to the critical value Pc, and T is 
an appropriate energy parameter. As a consequence the; 
entropy of each hot quasiparticle becomes anomalously 
large. This can be seen from the expression for entropy 
(5) per particle, 



N 



dl 
du) 



to tan 



Here /(w) is the Fermi function, and t{uj) is quasiparticle 
lifetime obtained from the inverse of imaginary part of 
self energy. From the above expression it is easy to see 
that S/N oc 1/Z. Over the hot region, keeping only the 
leading term, Z~^ ~ ln(l/J). Then, 



S/N oc Ar(o)r 



epiOs J y J 



(7) 



where M{0) is the density of states of the non-interacting 
system at the Fermi energy. For T > T{p~pc), the tem- 
perature dependence of entropy is S* oc nn(l/T), which 
is different from Fermi liquid behaviour {S oc T). This 
gives rise to the anomalous logarithmic temperature de- 
pendence of specific heat. In the past^'"^^ this behaviour 
has been understood from a purely bosonic point of view 
following the formalism of Hertz and Millis. For the 
spin fluctuations the Gaussian part of the action gives 
a free energy F oc T^lnT, which explains the ln(l/T) 
behaviour of C/T. Thus, here we find that there is agree- 
ment between the results of the spin-fermion model and 
the pure bosonic model. 



IV. THERMOPOWER 

From our discussion on entropy, it is natural to expect 
that this entropy enhancement should be seen in the mea- 
surement of thermopower (St). This is because one can 
think of thermopower as proportional to the correlation 
function between the heat current and the particle cur- 
rent, and heat current involves the transport of entropy 
due to temperature and electric potential gradients in the 
system. Strictly speaking, thermopower is defined as a 
ratio of two correlation functions, i.e. 



eTLn 



where 



Li2 = lim ^Im T dre^^ ( TJq (r) ^(0)), 

is the correlation function between heat current (jg) and 
particle current (j), and 



ill = 



lim — r-Im 



f 

Jo 



(ire--(T,j(r)-j(0)), 



is the correlation function between particle currents. Ln 
is a measure of electrical conductivity (a = e^Ln). Here 
we are ignoring the tensor nature of Ln and L12, and 
assuming that temperature and potential gradients and 
the thermal current are along the major symmetry di- 
rections of the lattice so that the tensors are diagonal. 
We express the single particle energies with respect to 
the chemical potential and assume that chemical poten- 
tial in the sample is uniform. The expression for heat 
current is given by. 



Jg 



-y 

p,cr 



Vp (( 



P,<t'-P>o- 



C^P,(tCp,(7^ • 



In principle, heat current will have a second term of 
the form | Y^^^^ Vfet/(k)(nk,<7«-k,(T - W-k.^^k.^), where 
C/(k) is the Fourier transform of the interaction term be- 
tween the electrons. However, such a term is quartic in 
fermionic operators and generates only subleading contri- 
butions in our calculation. We will also ignore corrections 
to the particle current and heat current vertices due to 
exchange of spin fluctuations. These vertex corrections 
are nonsingular, and change only the numerical prefac- 
tor (which we do not attempt to calculate) of our leading 
term, because the spin fluctuations are peaked around 
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a finite wave-vector. With these approximations the ex- 
pressions for the correlation functions can be re-expressed 
in a more transparent form as 



df 



Here Vp = dep/dp is the quasiparticle velocity, and 
A (p, uj) is the spectral function defined as 



Aip,Lu) 



(f-ep) +t{uj)- 



The evaluation of Ln is more straightforward and we 
will examine it first. The momentum sum can be con- 
verted into an integral over various energy surfaces. The 
dominant contribution is from the Fermi level, and we 
get 



We have already noted that over the hot region t(w) oc 



For the frequency integral since w T, we get 



Ln oc 



T 



(8) 



This result^^^ simply reiterates what we had noted be- 
fore, that when the hot carriers dominate transport. 
Act oc l/T. Now for L12, we first notice that the ex- 
pression is odd in frequency. This is because L12 is a 
measure of particle-hole asymmetry in the system. In 
our calculation we will consider as phenomenological in- 
put two different sources of such asymmetry. One such 
source is from the density of states, so that M{iS) = 
A/'(0)-hwA/''(0)-hO(wVe|), where TV' (0) 7^ only if there 
is particle-hole asymmetry in the bare non-interacting 
system of electrons. The second source of asymmetry will 
be from the quasiparticle lifetime which, for the hot car- 
riers, we write as t^^{lo) = {gQ/ep^^s) \^\ (1 + tuj). Here 
the second term is a possible particle-hole asymmetric 
term in scattering lifetime, r is a typical scattering time, 
and u) < r~^. After the energy integral around the Fermi 
surface we get. 



df 

duj(-^)ujT(ujW(Lu/Z) 
auj 



( EF^i 

I 9l 



vl {TN'{0)/Z + rA/'(0)T} . (9) 



The first term in the above equation is from the asym- 
metry in density of states, and the second term is from 
the asymmetry in quasiparticle lifetime. We note that 




FIG. 1: Magnetization (M) vs external magnetic field (H) 
at diflerent temperatures : a) T=0.15K, b) T=0.3K and c) 
T=0.8K. The discrete points are experimental. The solid lines 
are fits using equation (10). 



the factor of 1/Z, which leads to entropy enhancement, 
is associated with the asymmetry in density of states. 
Thus, the first term is the dominant one and eventu- 
ally gives anomalous temperature dependence to ther- 
mopower. For this leading term we can write 



(10) 



In the regime where T > r{p - Pc), St/T oc ln(l/T), as 
has been observed^*'^^ in thermopower measurement on 

CeCu6-a;AU:E- 



V. CONCLUSION 

To check the consistency of our model and calculation, 
we need to estimate the high energy scale (namely, Tj^) 
of CeCus.gAuo.i. For this purpose, we have fitted an ap- 
proximate form of the free energy function (F) that will 
match with the experimental results at low temperature 
and in the presence of magnetic field (H) . The function 
that matches well with the experiment has the form. 



F{T,H)/kB^X{T,H) In 



2 cosh 



\Y{T, H) 



(11) 
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FIG. 2: Entropy (S) per Ce atom vs temperature (T) at dif- 
ferent magnetic fields : a) H=OT, b) H=1.5T and c) H=3T. 
The discrete points are experimental, obtained by numerically 
integrating data from specific heat measurement. The solid 
lines are fits using equation (10). 

where 

X{T,H) - r*. + Ci 

Here Ci - C3 are parameters of the fitting function, fi is 
the effective magnetic moment of the Ce^+ ions in units 
of Bohr magneton (/is) and A = fis/kB — 0.67. We have 
chosen a simple possible form of the free energy which, at 
low temperatures (T <^ T^), is consistent with the criti- 
cal form of free energy that is suggested by renormaliza- 
tion group calculation for 2d spin fluctuations,^^ namely 
F oc T^\ii{To/T). At high temperatures (T > Tk) it 
matches smoothly to an impurity model where the 4/ 
cerium electrons act as Kondo impurities. The uniform 
magnetic susceptibility in this regime is Curie- Weiss like, 
with xiT) oc n'^/T. This temperature dependence is cut- 
off at T^, below which x ^ ? down to zero temper- 
ature. The fitting function is chosen such that at very low 
temperature (T ^ 0), x(T) -x(0) oc -T.^" This limiting 
behaviour agrees with the form x ~ oq + l/(ai -I- a2T) 
which Rosch et. al.^ used to fit susceptibility data up to 
1.4 K. We also find that susceptibility derived from equa- 
tion (10) can describe reasonably well (with a difference 



of at most twenty percent) the data'' up to 6 K. The varia- 
tion of entropy and magnetization as functions of temper- 
ature and magnetic field that one expects from the above 
free energy matches well with the experiments (Figures 
1,2). From the fit we estimate to be around 15 K, and 
/i ^ 2.6. In the absence of magnetic field the specific heat 
coefficient (7 = C/T) can be written as 7 = aln{TQ/T). 
From the fit we estimate a = 0.5 J/mol-K^ and Tq — 9.4 
K, which have comparable order of magnitudes with the 
experimentally measured values a — 0.6 J/mol-K^ and 
To — 5.3 K.^ The logarithmic behaviour in specific heat 
and thermopower in CeCu6-a;Aua; is observed around 1 
K, which is well below T^. The experimental fits and 
the estimates suggest that the spin-fermion model that 
we have been considering is consistent with the experi- 
mental data. 

We now discuss the limitations of our calculation. We 
have completely ignored the interaction between the spin 
fluctuations (the uq term). This is justified since this 
term is marginally irrelevant. In our calculation we con- 
sidered only the lowest order diagram in the perturbation 
series in terms of the spin-fermion coupling. However, we 
have examined the lowest order spin-fermion vertex cor- 
rection, and found that it is well-behaved close to the 
QCP (described in the appendix). So we believe that the 
qualitative features of our calculation will not be modi- 
fied by including higher order terms of the series. This 
is very different from what is found in the 2d-spin-2d- 
fermion model, where the spin-fermion vertex is singular 
indicating a potential breakdown of the approach. So, 
if the 2d-spin 3d-fermion model breaks down, there is no 
trace of this breakdown in perturbation theory. 

From our calculation we see that irrespective of 
whether the system is clean or dirty, if there is a large 
enough hot region in the system, then both specific heat 
and thermopower should show anomalous logarithmic 
temperature dependence. 

Since the microscopic origin of the 2d spin fluctuations 
is not known, our model seem to be a flne tuned one 
rather than one that is expected intuitively. It would be 
interesting to investigate the origin of the 2d magnetic 
coupling, and why most of the Fermi surface is hot by 
means of microscopic flrst principles calculations. This 
study should be supplemented by an investigation of the 
2d-3d dimensional crossover, to estimate the energy scale 
at which it is expected to occur. We notice that speciflc 
heat and resistivity measurements on YbRh2Si2^ seem to 
indicate that the model, with most of the Fermi surface 
hot, is quite valid for it. From this we can conclude that 
we expect to see the behaviour St/T oc ln(l/T) from 
thermopower measurement on YbRh2Si2, probably over 
a wider range of temperatures than the Ce-material. 
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APPENDIX A: SPIN-FERMION VERTEX 

Here we describe the calculation of the spin-fermion 
vertex and show that at the QCP (6 — > 0) the vertex is 
not singular. This is important because otherwise our 
perturbative calculation will break down at low temper- 
ature near the QCP. With a singular vertex, the coupling 
constant between the electrons and the spin fluctuations 



will get strongly rcnormalizcd at low energy. The quali- 
tative features of the theory will change, in particular the 
electron self energy. We will express the lowest order cor- 
rection to the bare spin-fermion coupling as g = go{l+r). 
Since we are interested only in the hot electrons and their 
low energy interaction with the spin fluctuations, we will 
calculate the vertex T with all external frequency zero. 
The expression for the vertex will then be, 

G(p2-hk,w)x(Q + ax-Fk,w). (Al) 

Here pi and p2 arc two hot points that are connected by 
wave- vector Q -|- ax. Expressing the linearized spectrum 
near the two hot points as eik and e2k, we can rewrite 
the above expression as 



r = 4g'oJ2 —7 \ • 

k -^0 ^ (7Q+k + '^^J (eik + e2k) (w + eik) (w + eak) 



It is easy to check by simple dimensional analysis that as cutoff in the momentum space. 
i5 ^ 0, the above expression is finite. As an estimate we 
get F cx g'^A^/^ / (e'p^wl^^), where A is a dimensionless 
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